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Abstract: We report on a non-perturbative study of two dimensional V = (2, 2) super QCD. 
Our lattice formulation retains a single exact supersymmetry at non-zero lattice spacing, and 
contains Nj fermions in the fundamental representation of a U{Nc) gauge group. The lattice 
action we employ contains an additional Fayet-Iliopoulos term which is also invariant under 
the exact lattice supersymmetry. This work constitutes the first numerical study of this 
theory which serves as a toy model for understanding some of the issues that are expected 
to arise in four dimensional super QCD. We present evidence that the exact supersymmetry 
breaks spontaneously when Nf < Nc in agreement with theoretical expectations. 
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1. Introduction 

In recent years a new approach to the problem of putting supersymmetric theories on the 
lattice has been developed based on discretization of a topologically twisted version of the 
continuum theory [1, 2, 3, 4, 5, 6]. ^ Initially the focus was on lattice actions that target 
pure super Yang-Mills theories in the continuum limit, in particular AA = 4 super Yang-Mills 
[12, 13, 14, 15, 16]. For alternative approaches to numerical studies of AA = 4 Yang-Mills see 
refs. [17, 18, 19, 20]. However in [21] [22] these formulations were extended to the case of 
theories incorporating fermions transforming in the fundamental representation of the gauge 
group and hence targeting super QCD. The starting point for these later lattice constructions 
is a continuum quiver theory containing fields that transform as bifundamentals under a 
product gauge group U{Nc) x U{Nf). After discretization these bifundamental helds connect 
two separate lattices and, in the limit that the U{Nf) gauge coupling is sent to zero, yield a 
super QCD theory with a global U{Nf) flavor symmetry. This construction is described in 
detail in section 3. The lattice action we have employed in this work includes an additional 
Fayet-Illopoulos term which, while invariant under the exact lattice supersymmetry, generates 
a potential for the scalar fields. It is straightforward to show that this yields a non-zero 
vacuum expectation value for the auxiliary field (D term supersymmetry breaking) if Nf < Ve¬ 
in section 4. we show the results from numerical simulations of this theory which support this 

'^The same lattice theories can be obtained using orbifold methods and indeed supersymmetric lattice 
actions for Yang-Mills theories were first constructed using this technique [7, 8, 9, 10] and the connection 
between twisting and orbifold methods forged in [11] 
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conclusion; we measure a non-zero vacuum energy and show that a light state - the Goldstino- 
appears in the spectrum of the theory if Nf < Nc- In contrast we show that vacuum energy 
is zero and this state is absent from the spectrum when Nf > Nc which is consistent with the 
prediction that the theory does not spontaneously break supersymmetry in that case. 

2. The starting point: twisted Q = 8 SYM in three dimensions 

We start from the continuum eight supercharge (Q = 8) theory in three dimensions which is 
written in terms of twisted fields which are completely antisymmetric tensors in spacetime 
under the twisted SO(3) group. The original two Dirac fermions reappear in the twisted 
theory as the components of a Kahler-Dirac field T = {ri,'ipaiXab,dabc) where the indices 
a,b,c = 1... 3. The bosonic sector of the twisted theory comprises a complexified gauge 
held Aa = Aa + iBa containing the original gauge held Aa and an additional vector held 
Ba- This additional held contains the three scalars expected of the eight supercharge theory 
which, being vectors under the R symmetry, transform as a vector held after twisting. The 
corresponding action S = S'exact + 5'ciosed where 

•Sexact — ~~2 QA = 4 Q / d^xTi Xabix)Bab{x) + vix) [Da, Da] + 4(x)d(x) , (2.1) 

S' S' y L ' ' 2 

•Sclosed = - ^ j d^xTr [9abcix)D[cXab]ix)] ■ (2.2) 

Here all helds are in the adjoint representation of a U{N) gauge group X = 
and we adopt an antihermitian basis for the generators Tq. Da and Da are the continuum 
covariant derivatives dehned in terms of the complexihed gauge helds as Da = da + Aa and 
Da = da + Aa- The action of the scalar supersymmetry on the helds is given by 

QAa = V’a 
QAa = 0 
Qi^a = 0 

QXab — D ab 

Qrj = d 

QBabc = 0 (2.3) 

Notice that we have included an auxiliary held d{x) that allows the algebra to be off-shell 
nilpotent = 0. This feature then guarantees that 5exact is supersymmetric. The equation 
of motion for this auxiliary held is then 

d{x) = [Da,Da] 


- 2 - 


(2.4) 



The Q-invariance of Sdosed follows from the Bianchi identity^ 


^abc'^cJ~ ab O’ 


(2.5) 


To discretize this theory we place all fields on the links of a lattice. This 3d lattice consists of 
the usual hypercubic vectors plus additional face and body links. In detail these assignments 
are 


continuum field 

lattice link 

Aa{x) 

X ^ X + d 

Aa{x) 

X -t- a —)• X 

ll^aix) 

X —>■ X -|- o 

Xab 

X + d + b ^ X 

r/(x) 

X —)• X 

d{x) 

X —)> X 

dabc 

X—I’X-l-a-l-ft-l-c 


The lattice gauge field will be denoted U^{x) in the following discussion. For the scalar 
fields d{x), rj{x) the link degenerates to a single site. Notice that the orientation of a given 
fermion link field is determined by the even/odd character of its corresponding continuum 
antisymmetric form. The link character of a field determines its transformation properties 
under lattice gauge transformations eg. Ua{x) —)■ G{x)Ua{x)G\x + a). To complete the con¬ 
struction of the lattice action it is necessary to replace continuum covariant derivatives by 
appropriate gauged lattice difference operators. The necessary prescription was described in 
[4], [21], [23]. It is essentially determined by the simultaneous requirements that the lattice 
difference agree with the continuum derivative as the lattice spacing is sent to zero and that 
it yields expressions that transform as the appropriate link field under lattice gauge transfor¬ 
mations. The lattice difference operators acting on a field where (±) corresponding to 

the orientation of the field^, are given by: 




bjx + a) - + h) 

Ua{x + ^Jx + a) - 


( 2 . 6 ) 

(2.7) 




( 2 . 8 ) 

(2.9) 


VK,. bAx) . 


f)l,f)2,...,6n ^ ^ ■’bl,b 2 ,...,b„ 

a Jbub 2 ,...,bAA = = - “)’ 


'T’)( ) ^(i) 
J) 


( 2 . 10 ) 

( 2 . 11 ) 


^Note that it is also possible to write the 3d action completely in terms of an Q-exact form without a 
Q-closed term by employing an additional auxiliary field Babe 

^Note that ipaix) and 9abc{x) originate from lattice site x and are, thus, positively oriented. Xab{x), however, 
terminates at lattice site x and this therefore assigned a negative orientation. 
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where b = ™ equations (2.6) to (2.9). For example the continuum derivative Daipb 

becomes 

= Ua{x)'ipbix + a) - 'ijjbix)Ua{x + b) ( 2 - 12 ) 

This prescription yields a set of link paths which, when contracted with the link field Xab{x), 
yields a closed loop whose trace is gauge invariant: 


Tr 


Xab{x) (Uaix)'tpb{x + a) - 'lpb{x)Uaix + b 


(2.13) 


It has the correct naive continuum limit provided that (in some suitable gauge) we can expand 
Ua{x) = In + Aa{x). The held strength on the lattice, Tab{x), is dehned using the forward 
difference operator as: 

Tab{x) = vA^Ub{x). (2.14) 

In lattice QCD the unit matrix arising in this expansion is automatic since the link helds take 
their values in the group. However the constraints of exact lattice supersymmetry require that 
the lattice gauge helds take their values, like the fermions, in the algebra. In this case the unit 
matrix can then be interpreted as arising from giving a vev to the trace mode of the original 
scalar helds Ba- This feature is required by lattice supersymmetry but is only possible because 
we are working with a complexihed U{N) theory - another indication of the tight connection 
between twisting and exact supersymmetry. It also implies that the path integral dehning 
the quantum theory will use a hat measure rather than the usual Haar measure employed 
in conventional lattice gauge theory. Such a prescription would usually break lattice gauge 
invariance but again complexihcation comes to the rescue since the Jacobian resulting from 
a gauge transformation of the DU measure cancels against an equivalent one coming from DU. 


We now show how to use this three dimensional lattice model to construct a two dimen¬ 
sional quiver theory while maintaining the exact lattice supersymmetry. 


3. Two dimensional quivers from three dimensional lattice Yang-Mills 

Consider a lattice whose extent in the 3-direction comprises just two 2d slices. Furthermore we 
shall assume free boundary conditions in the 3-direction so that these two slices are connected 
by just a single set of links in the 3-direction - those running from X 3 = 0 to X 3 = 1. Ignoring 
for the moment any fields that live on these latter links it is clear that the gauge group can 
be chosen independently on these two slices. We choose a group U (W) for the slice at X 3 = 0 
and U{Nf) at X 3 = 1 and will henceforth refer to them as the W and Nf lattices. Denoting 
directions on the 2d slices by Greek indices = 1,2 the fields living entirely on these lattices 
are given by 


Nc 

: 'i>{x) = 


— Inc + ^ 

(3.1) 

Nf 

: 'F(x) = 


, IjSJ J, ; d 

(3.2) 


- 4 - 




In these expressions x{x) denotes the coordinates on the Nc{Nf) lattice and denote 

the Nc{Nf) X Nc{Nf) unit matrix respectively. Now consider fields that live on the links 
between the and Nf lattice. These must necessarily transform as bi-fundamentals under 
U{Nc) X U{Nf). We have, 

Nc X Nj . ^bi-fund(3^) 3;) — (^■^) 

The second equality in the above equation is a mere change of variables and corresponds to 
labeling fields according to their two dimensional character. The complete field content of 
this model is summarized in the table below: 


W-lattice 

X 

Bi-fundamental fields 
(x,x) , (x,x) 

Nf-lattice 

X 

A^{x) 

(/>(x, x) 

Af,{x) 

r]{x) 

A(x, x) 

Vix) 


A^(x-k/x,x) 

^m(^) 

Xi^u{x) 

Xf_i,y{x,x + n + i^) 

Xiiiy{x) 


Defining G(x) as a group element belonging to U{Nc) and H(x) to U{Nf) the lattice gauge 
transformations for the bi-fundamental fields are as follows: 


4>{x) —> G{x)(j){x)H'^ {x) 

X{x) —> G{x)X{x)H^x) 

A^(x) H{x + fi)X^{x)G\x) 

Xfj.uix) G{x)Xfj,,y{x)H^x + + 1 ^) (3.4) 


It is crucial to note that this generalization of the original lattice super Yang-Mills theory to 
a quiver model is completely consistent with both the quiver gauge symmetries and the exact 
supersymmetry. For example the 3d term given in eqn. 2.13 yields a bi-fundamental term of 
the form 


Tr 


Xf,{x) (u^{x)X{x + fi) - X{x)Ufj_{x)^ 


(3.5) 


which is invariant under the the generalized gauge transformations given in eqn. 3.4. Thus, 
the above construction lends us a consistent lattice quiver gauge theory containing both 
adjoint and bi-fundamental fields transforming under a product U{Nc) x U{Nf) gauge group. 
Consider now setting the U{Nf) gauge coupling to zero. This sets = Inj up to gauge 
transformations and it is then consistent to set all other fields on the Nf lattice to zero. The 
original U{Nf) gauge symmetry now becomes a global U{Nf) flavor symmetry which acts on 
a set of complex scalar fields cj) transforming in the fundamental representation of the gauge 
group and their fermionic superpartners (A, A^, A^^.). The situation is depicted in figure 1. At 
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this point we have the freedom to add to the action one further supersymmetric and gauge 
invariant term - namely r^^Trd(x) = rQ^^Trry. This is a Fayet-Iliopoulos term. Its 
presence changes the equation of motion for the auxiliary field 

d{x) = + <i){x)^{x) - riN, (3.6) 

with I]\f^ a Nc X Nc unit matrix. The SUSY transformations for the remaining adjoint and 
fundamental fields are: 


Adjoint Fields 

Fundamental fields 


Q4> = A 


Q^ = 0 

QV'm = 0 

QA = 0 

Qx/j-i' ~ iJ,v 

QA^ = -Vfj^cp 

Qr] = d 

> 

II 

O 


After integration over d the Fayet-Iliopoulos term yields a scalar potential term which will play 
a crucial role in determining whether the system can undergo spontaneous supersymmetry 


(7(We) SYM Adjoint Model 



UiNp) SYM Adjoint Model 


Figure 1: 3d quiver model 

































breaking. The final action may be written as 


-^adj = -T^y{x)T^y{x) - - r]{x)vl ^ip^^{x) - 'iIj^]{x) 


(3.7) 


‘S'fund = ^ Tr -vy’(p{x)vj;^'>(j){x)-- {4>{x)(l){x) - li) + V^’Uf,{x) {(p{x)(j){x) - vi) 

X *- 

- \r]{x)X{x)'^{x) + |a^(x)pWA(x) - \^,{x)Tp^{x)ct){x + //)| 


| a ^ i .( x ) V^^'’\y{x) - \^y{x)4){x + // + i^)XfMuix)'^ 


(3.8) 


In practice we have also included the following soft SUSY breaking mass term, S'softj in the 
adjoint action, S'adj in equation (3.7): 


'S'soft A* 


-IV (U,U,) - 1 


(3.9) 


Such a term is necessary to create a potential for the trace mode of the twisted scalar fields 
as we have discussed earlier. In principle we should extrapolate —>• 0 at the end of the 
calculation and so we have obtained all our results for a range of In practice we observe 

that these soft breaking effects are rather small. 


Finally, the lattice coupling k appearing above is given by: 


^ NcLT 

2XA ■ 


(3.10) 


Here, A = g'^Nc is the dimensionful ‘t Hooft coupling, L and T are the numbers of points in 
each direction of the 2d lattice and H is a continuum area - the importance of interactions 
in the theory being controlled by the dimensionless combination XA. When we later discuss 
our numerical results we refer to this dimensionless combination as simply A. 


4. Vacuum Structure and SUSY Breaking Scenarios 

Let us return to the equation of motion for the auxiliary field d{x). If we sum the trace of 
this expression over all lattice sites and take its expectation value we find 

(^Tr d{x)) = [(t>{x)^{x) - vIn,)) (4.1) 

X X 

Since the lefthand side of this expression is the expectation value of the Q-variation of some 
operator the question of whether supersymmetry breaks spontaneously or not is determined 
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by whether the righthand side is non-zero. Indeed after we integrate over the auxiliary field 
d we find a scalar potential of the form 

Nf _ 2 

5'Dterm = ^Tr {x)4>\x) - , (4.2) 

Consider the case where Nf < Nc- Using SU{Nc) transformations one can diagonalize the 
Nc X Nc matrix cjxj). In general it will have Nf non-zero real, positive eigenvalues and Nc — Nf 
zero eigenvalues. This immediately implies that there is no conhguration of the fields cj) where 
the potential is zero. Indeed the minimum of the potential will have energy r‘^{Nc — Nf) and 
corresponds to a situation where Nf scalars develop vacuum expectation values breaking the 
gauge group to U{Nc — Nf). The situation when Nf > Nc is qualitatively different; now the 
rank of (jxp is at least Nc and a zero energy vacuum conhguration is possible. In such a sit¬ 
uation Nc scalars pick up vacuum expectation values and the gauge symmetry is completely 
broken. 

For the case when Nf < Nc where Q-supersymmetry is expected to break we would ex¬ 
pect the spectrum of the theory to contain a massless fermion - the goldstino [24]. To see 
how this works in the twisted theory consider the vacuum energy 

/ 0, (4.3) 

which is equivalent to < {Q,0} 0 for some operator O. In the two dimensional twisted 

theory the relevant part of the supersymmetry algebra is {Q, Q^} = [25] so that eqn. 4.3 

is equivalent to 

(0|{Q,Qo}|0)/O, (4.4) 

Note that the equation above involves both the scalar Q and the I-form supercharge Q^. 
Corresponding to these supercharges are a set of super cur rents, J and whose form can 
be derived in the usual manner by varying the continuum twisted action under infinitesimal 
spacetime dependent susy transformations. This yields gauge invariant supercurrents on the 
lattice of the following form 

Jix) = d{x) + ..., (4.5) 

Jo{x) = r]{x)d{x) + ..., (4.6) 

and using the equations of motion, the auxiliary field d(x) can be replaced by 

d{x) = ^ + [(l}{x)'(j)ix) -tInc] (4.7) 

fi=l,2 

We therefore expect a possible Goldstino signal to manifest itself in the contribution of a light 
state to the two-point function: 


Cit) = {0\0{x)0\y)\0) 


( 4 . 8 ) 


where ‘t’ corresponds to [x^ — y^) and a suitable set of lattice interpolating operators are 
given by: 


0{x) = Tr 


'^'il^f,{x)U^{x) {(j){x)(t){x) 


riNc) 


(4.9) 


L 


and 

0'{y) = Tr [r]{y) {(t){y)^{y) - rlN,)] • (4.10) 


5. Numerical Results 

We employ a RHMC algorithm to simulate our system having first replaced all the twisted 
fermions in our model by corresponding pseudofermions - see for example [26] [27]. The 
simulations are performed by imposing anti-periodic (thermal) boundary conditions on the 
fermions along one of the two space-time directions. This is done to avoid running into the 
fermion zero modes resulting from the scalar component of the twisted fermion, rj. As dis¬ 
cussed in [28] [29] this has the added benefit of ameliorating the sign problem for these lattice 
theories. This breaks supersymmetry explicitly by a term that vanishes as the lattice volume 
is increased. 

In this section, we contrast results from simulations with Nf = 2, Nc = 3 corresponding 
to the predicted susy breaking scenario with results from simulations with Nf = 3, Nc = 2 
- the susy preserving case. We ran our simulations for three different values of the ‘t Hooft 
coupling, A = 0.5,1.0 and 1.5 and observed the same qualitative behavior for the different 
values of A. The results presented in this section correspond to A = 1.0. The FI parameter, 
r, is a free parameter and is set to 1.0 for the rest of the discussion. 

As a first check, we compared the expectation value of the bosonic action with the theo¬ 
retical value obtained using a supersymmetric Ward identity 

< AtSboson >= + NcNf^ V. (5.1) 

In appendix A. we show how to compute this value. Figure 2 shows a plot of the bosonic 
action for various values of the soft SUSY breaking coupling [i. In principle we should take 
the limit /r —?■ 0 although it should be clear from the plot that the /r dependence is in fact 
rather weak. We have normalised the data to its value obtained by assuming supersymmetry 
is unbroken. The red points at the bottom of the figure denote the SUSY preserving case and 
it can be observed that they agree with the theoretical prediction. This is to be contrasted 
with the case when Nf < Nc denoted by the blue points which shows a large deviation from 
eqn. 5.1 and is the first sign that supersymmetry is spontaneously broken in this case. 

The spatial Polyakov lines shown in figure 3 also show a distinct difference between the 
Nf < Nc and Nf > Nc cases. The red lines where |P| ~ 1 correspond to the SUSY preserving 
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Figure 2: Normalized bosonic action vs soft breaking coupling /r for A = 1.0 for a 16x6 lattice 


case and are consistent with a deconfined or fully Higgsed phase. Indeed the Polyakov line is a 
topological operator and in a susy preserving phase should be coupling constant independent 
consistent with what is seen. The blue line in the lower half of the plot corresponds to smaller 
values which is qualitatively consistent with the predicted partial Higgsing of the gauge field 
in the phase where supersymmetry is spontaneously broken. 

One of clearest signals of supersymmetry breaking can be obtained if one considers the equa¬ 
tion of motion for the auxiliary field eqn. 4.2. We expect the susy preserving case to obey 

[<^(x)^(x)] = 1. (5.2) 

The red points, corresponding to {Nf > Nc) are consistent with this over a wide range of 
II. We attribute the small residual devaition as /i —)• 0 to our use of antiperiodic boundary 
conditions which inject explicit Q susy breaking into the system. The simulations with 
Nf < Nc (blue points) however show a clear signal for spontaneous supersymmetry breaking 
with the value of this quantity deviating dramatically from its supersymmetric value even as 
/r —0. 

Finally we turn to our results for a would be Goldstino. We search for this by computing the 
following two point correlation function 

C{t) = '^< 0'{y, t)0{x, 0) > (5.3) 

^,y 
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Figure 3: Spatial Polyakov line vs p, for A = 1.0 on an 16x6 lattice 
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Figure 4: ^Tr00 vs for a’t Hooft coupling of A = 1.0 on an 16x6 lattice 


where 0^{y^t) and O(x,0) are fermionic operators given by: 

0{x, 0) = ip/j.ix, 0)Uij,{x, 0) [(j){x, 0)4>{x, 0) - tIatJ (5.4) 

0'{y, t) = r]{y, t) [(l){y, t)^{y, t) - tIatJ . (5.5) 
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A = 1.0 ;/i = 0.3 



Figure 5: Correlation function C(t) for A 


1.0 and /i = 0.3 on various asymmetric lattices 


Since it is computationally very cumbersome to evaluate the above correlation function for 
every lattice site x at the source we instead evaluate the correlator for every lattice site y for 
a few randomly chosen source points x. In figure 5 we show the logarithm of this correlator 
as a function of temporal distance for a range of spatial lattice size, L = 6,8,12 and 14. 
The anti-periodic boundary condition is applied along the temporal direction corresponding 
to T=16 and for both Nf > Nc and Nf < Nc- The approximate linearity of these curves is 
consistent with the correlator being dominated by a single state in both cases. However when 
Nf > Nc the amplitude of this correlator is strongly suppressed relative to the case where 
Nf < Nc- Furthermore the effective mass extracted from fits to this latter correlator (figure 
6 ) falls as the spatial lattice size (L) increases, consistent with a vanishing mass in the large 
volume limit. The lines in figure 6 show fits to 1/L - the smallest mass consistent with the 
boundary conditions - the dashed green line is a ht constrained to go through the origin while 
the dotted red line allows the intercept to float. This is just what we would expect of a would 
be Goldstino arising from spontaneous breaking of the exact Q-symmetry. 

6. Conclusions 

In this paper, we have reported on a numerical study of super QCD in two dimensions. The 
model in question possesses M = (2, 2) supersymmetry in the continuum limit while our 
lattice formulation preserves a single exact supercharge for non zero lattice spacing. It is 
expected that the single supersymmetry will be sufficient to ensure that full supersymmetry 
is regained without fine tuning in the continuum limit. This constitutes the first lattice study 
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Figure 6: Goldstino mass derived from fits Meff vs inverse transverse lattice size, L ^ 


of a supersymmetric theory containing fields which transform in both the fundamental and 
adjoint representations of the gauge group. Our lattice action also contains a Q-exact Fayet- 
Iliopoulos term which yields a potential for the scalar fields. The lattice theory possesses 
several exact symmetries; U{Nc) gauge invariance, Q-supersymmetry and a global U{Nf) 
flavor symmetry. 

It is expected that the system will spontaneously break supersymmetry if Nj < Nc- The 
arguments that lead to this conclusion depend on the inclusion of the Fayet-Iliopoulos term. 
Such a term is rather natural in our lattice model since the formulation requires U{Nc) gauge 
symmetry. Notice, though, that the free energy of the lattice model does not naively depend 
on the coupling r as long as it is positive since the Fayet-Iliopoulos term is Q-exact.^ Our 
numerical work is fully consistent with this picture; we have examined several supersymmetric 
Ward identities which clearly distinguish between the Nf < Nc and Nf > Nc situations and 
we have observed a would be Goldstino state in the former case. 

There are many directions for future work; inclusion of anti-fundamentals fields is straight¬ 
forward since it merely corresponds to including the bifundamental fields truncated from the 
A^j-lattice. Observations of phase transitions in such models as the parameters are varied can 
then potentially probe sigma models based on Calabi-Yau hypersurfaces [30]. It is possible 
that the SU{N) theories could be studied by deforming the moduli space of the lattice theory 

contrast for r < 0 we would expect supersymmetry breaking for any value of Nf/Nc- Thus one expects 
a phase transition in the Nf > Nc theory at r = 0. 


- 13 - 









using ideas similar to those presented in [31]. This would allow direct contact to be made to 
the continuum calculations of Hori and Tong [32], Finally the lattice constructions discussed 
in this paper generalize [33] to three dimensional quiver theories leaving open the possibility 
of studying 3D super QCD using lattice simulations. 


A. Calculating the Bosonic Action 

Consider the partition function 

Z = [ IlXe"'"(2A+Sc) 


(A.l) 


where DX denotes the measure over all boson and fermion fields and Sc the Q-closed term. 
We start by rescaling the field Oabc i^(^abc to remove the coupling n from in front of the 
Q-closed term. This yields 




(A.2) 


with V the two dimensional volume. Notice that NcNf is the number of fermions at each site 
resulting from the 3d 9 held. Differentiating with respect to n gives 


dlnZ _ NcNfV dlnZ' 

dn K Sk 


(A.3) 


The last term in the righthand side being Q-exact would yield zero in the original theory 
containing a d-held. However in the action we simulate this held is integrated out yielding 
instead a contribution of — Putting these pieces together we hnd 


K< Sb> +K< Sf>= -NcNfV - 


(A.4) 


The expectation value of the fermionic action can be gotten by scaling arguments since the 
fermions occur only quadratically in the action yielding 


4H 


K<Sf>=-—{N^ + NcNf) 
Collecting terms yields the hnal result quoted previously 


K<Sb>=V Qa2 + NcNf^ 


(A.5) 


(A.6) 
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